mi  null 


saeNce  ApplicauoNs,  mc. 


gi  s 


ol  174 


■  U’tf  njp. r  Ckj*.** 


SAI-131-81-275-LJ 


JIJ 


ANNUAL  REPORT 
CONTRACT  N00014-77-C-0489 
SATELLITE  MEASUREMENTS  OF  ATMOSPHERIC  AEROSOLS 


August  26,  1981 


Principal  Investigator: 


M.  Griggs 

Science  Applications,  Inc. 
Electronic  Vision  &  Systems  Division 
11526  Sorrento  Valley  Road,  Suite  A 
San  Diego,  California  92121 


Prepared  for: 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  Virginia  22217 


_ _ STS?”*3" 


SCIENCE  APPLICATIONS.  LA  JOLLA,  CALIFORNIA 

ALBUQUERQUE  •  ANN  AR80R  •  ARLINGTON  •  ATLANTA  •  BOSTON  •  CHICAGO  •  HUNTSVILLE 
LOS  ANGELES  •  McLEAN  •  PALO  ALTO  •  SANTA  BARBARA*  SUNNYVALE  •  TUCSON 


P.O.  Box  2361,  1200  Pmpict  Strait,  La  Jolla.  California  92037 


SECURITY  CLASSIFICATION  of  This  PACE  Qua  Enfred) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  iPEFO RT  WtrWBEO  -  2.  GOVT  ACCESSION  NO. 

'  / SAI-131  -81-275-LJ  '  Ah-A/^  4 ^ 

>  RECIPIENT'S  CATALOG  NUMBER 

y  ’■  <  * 

4.  TITL.  C  (end  Subtitle) 

Satellite  Measurements  of  Atmospheric  Aerosols  ^ 

[  Annual  Rep«rf*lHovember  15^ 
TOSITto  August  26,  1981 

C.  PERFORMING  ORG.  REPORT  NUMBER 

SAI-131-81-275-LJ 

T.  Au  T  HORf  aj 

M.  Griggs 

*.  CONTRACT  OR  GRANT  NUMBCRf*; 

'^N00014-77-C-0489  y" 

9.  PERFORMING  ORGANI2ATION  N»ME  AND  ADDRESS 

Science  Applications,  Inc. 

P.  0.  Box  2351 

La  Jolla,  CA  92037 

10.  PROGRAM  ELEMENT.  PROJECT.  TACK 
AREA  &  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research,  Code  465  / 

800  North  Quincy  Street  •• 

Arlington,  VA  22217 

1*.  REPORT  OATE 

.  Augwafc  26* JU81 

t*.  NUMBER  OF  PAGES 

_ is _ /~ }  :  _ 

14.  MONITORING  AGENCY  NAME  6  ADDRESSflf  dIHrront  Irom  Controlling  Ollier) 

15.  SECURITY  CLASS,  (of  thte  report) 

UNCLASSIFIED 

IS*.  OECL  ASSlFlC  AT'CN  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ol  1  hit  Report; 

Approved  for  public  release;  distribution  unlimited. 

*7.  DISTRIBUTION  STATEMENT  ( © /  the  mbetreet  entered  tn  Block  20,  if  different  from  Report) 

Same  as  report. 

ie.  supplementary  notes 

is  KEY  WORDS  (Continue  on  reveree  eide  if  neceesmry  end  Identify  by  block  number) 

Maritime  Aerosols  Aerosol  Scattering 

Satellite  Measurements  Sun  Photometer  Measurements 

Error  Analysis 

Ocean  Radiance 

>4^.  ABSTRACT  (Confirm*  on  reveree  eide  ff  neceeeery  end  Identify  by  block  number) 

Y Analysis  of  N0AA-6  AVHRR  data  and  ground- truth  measurements  at  two  locations  has 
found  a  radiance-aerosol  content  relationship  which  agrees  very  well  with  that 
previously  determined  for  Landsat  2  data.  A  new  technique  was  developed  to  use 
AVHRR  Channels  1  and  2  radiances  to  infer  aerosol  size  distribution  information, 
in  addition  to  the  optical  thickness.  The  technique  was  successfully  tested 
with  a  small  data  set  obtained  in  a  ground-truth  experiment  on  the  USNS  Hayes. 

A- 

\ 

DD  1473  EOITION  OF  I  NOV  tS  IS  OBSOLETE 

i 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  f»>i»n  DM  Entrrrd) 


ABSTRACT 


Analysis  of  NOAA-6  AVHRR  data  and  ground-truth  measurements  at  two 
locations  has  found  a  radiance-aerosol  content  relationship  which  agrees 
very  well  with  that  previously  determined  for  Landsat  2  data.  Problems  were 
again  found  in  the  radiometric  calibration  of  the  satellite  sensor,  but  com¬ 
parison  of  Channels  1  and  2  radiance  values  resulted  in  correction  factors 
being  determined  for  both  channels.  A  new  technique  was  developed  to  use 
AVHRR  Channels  1  and  2  radiances  to  infer  aerosol  size  distribution  informa¬ 
tion,  in  addition  to  the  optical  thickness.  The  technique  was  successfully 
tested  with  a  small  data  set  obtained  in  a  ground-truth  experiment  on  the 
USNS  Hayes. 
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1.  INTRODUCTION 


(1  2) 

The  investigation  of  the  satellite  technique'  '  to  measure 
tropospheric  aerosols  over  the  ocean  has  continued  both  experimentally  and 
theoretically.  Ground- truth  measurements  and  NOAA-6  satellite  measurements 
have  been  obtained  at  San  Diego  and  on  the  USNS  Hayes  in  the  Atlantic.  These 
data  have  been  analyzed  and  compared  with  the  earlier  Landsat  2  results. 
Theoretical  calculations  were  made,  and  led  to  the  preliminary  development 
of  a  new  method  to  use  the  NOAA-6  AVHRR  Channels  1  and  2  radiance  data  to 
infer  information  on  the  aerosol  size  distribution  in  addition  to  the  optical 
thickness.  This  has  great  potential  in  the  use  of  satellite  data  to  routinely 
predict  the  useful  range  for  FUR  systems,  and  other  electro-optical  systems, 
over  the  oceans. 

The  further  development  of  this  technique  can  ultimately  lead  to  a 
computer  code  to  automatically  convert  the  measured  satellite  radiances  to 
global  ocean  maps  of  aerosol  optical  thickness  (as  a  function  of  wavelength), 
surface  visibility,  predicted  slant  range,  and  a  size  distribution  parameter 
if  desired.  In  addition  to  the  routine  production  of  daily  global  aerosol 
maps  for  Naval  operations,  the  data  could  be  used  to  build  a  climatology  base 
for  aerosol  properties  over  the  oceans  for  use  in  Naval  systems  planning. 
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2.  SATELLITE  AND  GROUND-TRUTH  MEASUREMENTS 


Ground-truth  data  were  acquired  at  the  time  of  the  NOAA-6  satellite 
overpasses  (approximately  0730  l.s.t.),  at  San  Diego  and  on  the  USNS  Hayes 
off  the  Virginia  coast.  The  Advanced  Very  High  Resolution  Radiometer  (AVHRR) 
on  the  NOAA-6  has  four  channels,  two  of  which  are  useful  for  this  investiga¬ 
tion,  and  provides  daily  coverage  with  1  km  spatial  resolution.  The  useful 
channels  for  this  study  are  Channel  1,  at  0.65  urn  (almost  identical  to  the 
Landsat  MSS5),  and  Channel  2,  at  0.85  pm. 

Four  sets  of  useful  satellite-ground  truth  coincidences  were  ob¬ 
tained  on  the  USNS  Hayes,  and  seventeen  sets  were  obtained  at  San  Diego.  The 
sunphotometer  observations  on  the  USNS  Hayes  included  measurements  at  0.5  pm 
and  0.88  pm,  thus  permitting  a  direct  comparison  of  the  size  distribution 
parameter  determined  from  the  AVHRR  radiances  at  0.65  pm  and  0.85  pm.  The 
San  Diego  ground- truth  unfortunately  had  useful  measurements  only  at  0.5  pm, 
so  that  the  satellite  values  of  size  distribution  could  not  be  checked. 


3.  THEORETICAL  CALCULATIONS 


(2) 

It  was  shown  previously'’  '  that  the  radiance-aerosol  content  re¬ 
lationship  varies  with  the  size  distribution,  which  is  typically  well  repre¬ 
sented  in  the  real  atmosphere  by  the  Junge  distribution: 

dn(r)  =  Cr"v  d  log  r  (cm""*)  (1) 

where  n(r)  is  the  number  of  particles  with  radius  r,  and  C  is  a  constant 
depending  on  the  number  of  particles  per  unit  volume. 

The  previous  calculations  showed  that  the  spectral  variation  in 
the  visible-near  infrared  region  of  the  upwelling  radiance  in  nadir  viewing 
is  sensitive  to  the  aerosol  size  distribution.  The  earlier  Landsat  data 
were  not  sensitive  enough  to  investigate  the  use  of  satellite  spectral  mea¬ 
surements  to  infer  the  size  distribution,  but  the  AVHRR  on  NOAA-6  has  two 
channels  in  the  visible-near  IR  region,  and  is  a  10-bit  system,  which  should 
be  sensitive  to  the  small  radiance  changes  due  to  changes  in  the  size  dis¬ 
tribution.  Thus  calculations  were  made  so  that  a  two-channel  analysis  of  the 
AVHRR  radiance  data  could  be  made  in  order  to  infer  v  in  addition  to  the 
aerosol  content  (N). 

Calculations  were  made  for  v  =  2.5  and  v  =  3.5  at  0.85  pm  for 
aerosol  contents  N  and  2N  at  just  four  sun  zenith  angles  (60°,  66°,  72°  and 
78°).  Similar  calculations  were  made  at  0.65  pm  for  v  =  2.5  to  supplement 
the  existing  detailed  0.65  pm  calculations  for  v  =  3.5. 

In  order  to  infer  v  from  the  satellite  data  the  Channels  1  and  2 
radiances  (after  correcting  the  Channel  2  radiance  for  water  vapor  absorp¬ 
tion)  are  compared  with  the  theoretical  radiances  for  the  sun  and  view  angles 
at  the  time  of  the  measurement.  Model  values  of  N  and  v  are  chosen  so  that 
the  model  radiances  agree  with  the  measured  radiances  in  both  AVHRR  Channels. 
Thus  to  analyze  the  satellite  data,  interpolation/extrapolation  must  be  made 
for  the  non-tabulated  values  of  sun  angle  and  v.  This  is  currently  done 


manually,  and  assumes  a  linear  variation  of  the  radiances  with  both  parameters. 
This  procedure  must  be  computerized  if  more  extensive  satellite  data  are  ana¬ 
lyzed  in  the  future. 

The  ground-truth  measurements  of  the  aerosol  optical  thickness  (r^) 
may  be  used  to  determine  a  value  of  \>  for  comparison  with  the  satellite  value 
if  multispectral  sunphotometer  data  are  available.  If  a  Junge  distribution 
is  assumed,  it  has  been  shown  (e.g.,  Bullrich^),  to  a  close  approximation, 
that: 

ta  =  BA-(v“2>  (2) 


where  B  is  a  constant. 

Thus  values  of  at  two  wavelengths  can  be  used  to  calculate  v. 
These  wavelengths  should  be  the  same  as  those  of  the  satellite  measurements 
for  an  ideal  ground- truth  experiment. 


4.  DATA  ANALYSIS 


The  AVHRR  data  obtained  at  San  Diego  and  on  the  USNS  Hayes  have 
been  analyzed  using  the  table-look-up  code  developed  previously^  for  inter¬ 
preting  just  the  Channel  1  radiance,  and  by  the  two-channel  analysis  discussed 
in  Section  3. 

4.1  CHANNEL  1  ANALYSIS 

The  results  for  San  Diego  and  the  USNS  Hayes  are  given  in  Fig.  1 
which  shows  that  the  measured  radiances  (normalized  to  nadir  viewing  with  a 
sun  zenith  angle  of  63°)  are  lower  than  anticipated,  based  on  the  MSS5  regres¬ 
sion  line  for  San  Diego.  However,  it  is  believed  that  these  radiances  are 
low  due  to  radiometric  calibration  errors  in  the  AVHRR  sensor,  based  on  a 
comparison  of  the  radiances  measured  in  Channels  1  and  2.  It  is  recalled^’2^ 
that  the  radiometric  calibration  of  sensors  has  been  a  problem  in  intercompar¬ 
ing  the  radiance-aerosol  content  relationship  for  different  satellites. 

Figure  2  shows  the  Channel  1  radiances  plotted  against  the  Channel 
2  radiances  for  the  same  ocean  pixel  without  any  correction  for  sun  or  view¬ 
ing  angles  (the  Channel  2  radiance  is  corrected  for  water  vapor  absorption 
as  discussed  below).  It  is  seen  that  the  mean  curve  does  not  pass  through 
the  origin,  suggesting  that  the  AVHRR  radiance  calibrations  are  in  error, 
since  both  channels  should  indicate  zero  radiance  at  the  same  time.  Theoreti¬ 
cal  radiance  values  are  also  plotted  in  Fig.  2  for  the  same  sun  and  viewing 
angles  of  the  measured  radiances,  and  indeed  indicate  that,  regardless  of  the 
aerosol  content  and  size  distribution,  the  mean  curves  can  be  extended  through 
the  origin.  The  middle  of  the  three  dashed  curves  in  Fig.  2  show  where 
radiances  for  -0.5N  should  theoretically  lie.  The  mean  of  the  measured  data 

(the  solid  curve)  can  be  shifted  to  this  location  by  adding  0.5  mw/cm  /um/sr 

2 

to  Channel  1  radiances,  and  by  subtracting  0.25  mw/cm  /ym/sr  from  Channel  2 
radiances.  It  is  important  to  note  that  these  corrections  were  derived  inde¬ 
pendent  of  any  radiance-aerosol  content  relationship. 
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Aerosol  Content  (Ground  Truth) 


Figure  2.  Theoretical  and  Measured  AVHRR  Channel  1  Radiances 
Versus  Channel  2  Radiances. 
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With  the  corrected  values  of  the  Channel  1  radiances,  it  is  seen 
in  Fig.  1  that  normalized  radiances  show  excellent  agreement  with  the  Land- 
sat  relationship.  It  is  noted  that  a  group  of  five  data  points  with  N<0.35 
shows  radiances  greater  than  expected.  These  measurements  were  obtained  for 
sun  zenith  angles  greater  than  73°  where  the  flat  earth  model  used  in  deve¬ 
loping  the  table-look-up  code  was  expected  to  introduce  uncertainties. 

4.2  TWO  CHANNEL  ANALYSIS 

A  preliminary  investigation  has  been  made  into  the  use  of  AVHRR 
Channels  1  and  2  to  infer  N  and  v,  as  discussed  in  Section  3,  for  the  San 
Diego  and  USNS  Hayes  data.  The  results  of  this  analysis  are  given  in  Figs. 

3  and  4  and  in  Table  1. 

Results  for  only  11  San  Diego  data  sets  are  shown  since  the  results 
for  the  sun  zenith  angles  greater  than  73°  are  considered  unreliable  as  ex¬ 
plained  in  Section  4.1.  (Channel  2  satellite  data  was  also  unavailable  for 
one  of  the  other  overpasses.) 

The  results  in  Fig.  3  shows  good  agreement  between  the  N  values 
determined  from  the  satellite  and  the  ground- truth,  with  indication  that  the 
satellite  may  slightly  underestimate  the  larger  N  values.  The  reasons  for 
this  could  lie  in  the  assumptions  made  in  the  table-look-up  code,  in  the  sun- 
photometer  accuracy,  or  in  the  satellite  radiance  accuracy;  with  this  limited 
data  set  it  is  not  possible  to  presently  identify  the  reason. 

The  results  In  Fig.  4  show  excellent  agreement  between  the  v  values 
determined  from  the  satellite  and  the  ground-truth  on  the  USNS  Hayes.  The 
two  channel  analysis  was  also  carried  out  on  the  San  Diego  data  to  infer 
satellite  values  of  N  and  v,  but  unfortunately  no  multi  spectral  ground-truth 
data  were  available  to  compute  v  values  for  comparison.  However,  as  shown  in 
Table  1,  the  satellite  values  of  v  generally  range  between  3.5  and  4.7  which 
are  values  typically  measured  by  other  methods;  even  the  two  high  values  of 
6.0  and  6.5  have  been  measured  (by  other  techniques)  occasionally  by  other 
workers. it  is  noted  that  the  two  channel  analysis  results  in  values  of  N 
different  from  those  of  the  Channel  1  analysis,  but  overall  does  not  seem  to 
improve  on  the  agreement  with  the  ground-truth  values. 


•  San  Diego 
x  USNS  Hayes 


Ground  Truth 


Figure  3.  Comparison  of  AVHRR  and  Ground-Truth 
Measurements  of  Aerosol  Content. 


Figure  4.  Comparison  of  AVHRR  and  Ground  Truth  Measurements  of 
Aerosol  Size  Distribution.  (USNS  Hayes) 
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Table  1.  Comparison  of  Channel  1  Analysis  and 
Two  Channel  Analysis  at  San  Diego. 


Date 

N 

(Ground  Truth) 

2  Channel 

Channel  1 
(Corrected)  Only 

NSat 

vSat 

NSat 

6-13 

0.59 

0.66 

4.0 

0.47 

6-14 

0.85 

0.73 

3.5 

0.68 

7-11 

0.84 

0.65 

3.8 

0.67 

7-15 

0.56 

0.50 

6.0 

0.47 

7-16 

0.64 

0.56 

3.7 

0.57 

7-29 

1.37 

1.13 

4.7 

1.29 

7-30 

1.40 

1.10 

3.6 

1.09 

9-16 

0.86 

0.64 

4.2 

0.67 

9-17 

0.50 

0.54 

4.7 

0.50 

9-18 

0.55 

0.46 

2.3 

0.57 

10-18 

0.53 

0.70 

5.0 

0.98 

10-19 

0.27 

0.46 

6.5 

0.60) 

>0  >  73° 

10-20 

0.27 

0.64 

4.8 

0.507  0 
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4.2.1  Mater  Vapor  Correction  for  Channel  2 

The  spectral  bandpass  of  AVHRR  extends  from  0.67  um  to  1.2  pm 
with  50%  response  at  0.71  pm  and  1.02  pm.  This  spectral  region  includes 
several  water  vapor  absorption  bands^,  so  that  the  radiance  observed  by 
the  AVHRR  Channel  2  is  influenced  by  the  water  vapor  in  the  atmosphere. 
This  problem  has  previously  been  addressed  for  Landsat  data  since  the  MSS7 
bandpass  covers  part  of  the  same  spectral  region.  Pitts  et  al^  calcula¬ 
ted  the  atmospheric  transmission  for  the  MSS7  bandpass  as  a  function  of 
water  vapor  content,  based  on  high  spectral  resolution  calculations. 
Griggs'  '  showed  that  these  results  could  be  fit  to  a  good  approximation 
by  the  equation: 


t  =  exp(-.094u‘43)  (3) 

where  t  is  the  path  transmission  and 

u  =  w(l  +  sec  0q)  (4) 

where  w  is  the  vertical  water  vapor  content  in  precipitable  cm  and  eQ  is 
the  sun  zenith  angle. 

The  water  vapor  correction  was  applied  to  the  MSS7  radiance  values 
by  multiplying  them  by  1/t.  This  approach,  used  by  Pitts  et  al  is  a  simpli¬ 
fied  one,  ignoring  scattering  effects,  but  is  justified  since  the  correction 
is  small.  The  Landsat  study^  showed  that  the  transmission  factor  to  be 
applied  to  the  MSS7  radiances,  in  the  absence  of  water  vapor  information,  is 
0.81  with  an  uncertainty  of  +0.1. 

For  AVHRR  Channel  2,  using  the  spectral  information  given  by  Koepke 
P  and  Quenzel^,  the  transmission  factor  is  found  to  be  0.86  ±  0.07,  with  the 

actual  variation  of  transmission  being  approximated  by: 

t  =  exp(-.069u*43)  (5) 
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The  transmission  factor  in  the  AVHRR  bandpass  is  higher  than  in 
MSS7,  even  though  it  includes  an  additional  water  vapor  band,  since  the 
water  vapor  bands  represent  a  smaller  fraction  of  the  broader  bandpass. 

In  the  results  presented  above  in  this  section,  the  Channel  2 
radiances  for  San  Diego  were  corrected  using  a  transmission  factor  of  0.86, 
since  water  vapor  data  were  not  available.  For  the  USNS  Hayes  data  analysis, 
preci pi  table  water  vapor  maps  of  the  United  States  were  obtained  from  the 
National  Climatic  Center  in  Ashville,  and  used  to  determine  a  value  of  w. 

For  the  four  USNS  Hayes  data  sets,  w  varied  between  0.94  and  1.24,  result¬ 
ing  in  t  values  of  0.86,  0.87,  0.87  and  0.89,  which  are  close  to  the  fixed 
value  of  0.86  used  for  the  San  Diego  data. 

Thus,  it  is  believed  that  a  water  vapor  correction  factor  of  0.86 
can  be  satisfactorily  used  in  future  analysis  of  AVHRR  Channel  2  radiances. 
Based  on  the  calculations  on  the  radiance  variation  with  v  reported  prev¬ 
iously^,  the  uncertainty  of  ±0.07  in  the  water  vapor  correction  factor 
results  in  an  uncertainty  of  about  ±0.15  in  v.  A  more  detailed  error  anal¬ 
ysis  of  the  two  channel  analysis  is  required  to  determine  the  relative 
significance  of  this  uncertainty  in  comparison  with  those  due  to  errors  in 
the  satellite  data. 


5.  CONCLUSIONS 


Analysis  of  the  NOAA-6  AVHRR  data  and  the  ground-truth  measure¬ 
ments  at  two  locations  showed  that  the  radiance-aerosol  relationship  agrees 
very  well  with  that  previously  found  for  Landsat  2.  The  analysis  again  re¬ 
vealed  a  problem  in  the  calibration  of  the  satellite  sensor  as  reported 
previously  for  the  Landsat  and  GOES  satellites. 

A  new  technique  has  been  developed  to  infer  aerosol  size  distri¬ 
bution  information,  in  addition  to  the  aerosol  optical  thickness,  from 
analysis  of  the  AVHRR  Channels  1  and  2  radiances.  Application  of  this 
method  to  the  limited  data  set  from  the  USNS  Hayes  in  the  Atlantic  showed 
excellent  agreement  between  the  satellite  and  ground- truth  values  of  the 
size  distribution  parameter.  Further  ground-truth  data  are  required  to 
validate  the  technique,  and  calculations  are  needed  to  perform  an  error 
analysis. 
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